Are the renormalized band widths in TTF-TCNQ of structural or electronic origin? 

— An angular dependent NEXAFS study. 
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We have performed angle- dependent near-edge x-ray absorption fine structure measurements in 
the Auger electron yield mode on the correlated quasi-one-dimensional organic conductor TTF- 
TCNQ in order to determine the orientation of the molecules in the topmost surface layer. We find 
that the tilt angles of the molecules with respect to the one- dimensional axis are essentially the same 
as in the bulk. Thus we can rule out surface relaxation as the origin of the renormalized band widths 
which were inferred from the analysis of photoemission data within the one-dimensional Hubbard 
model. Thereby recent theoretical results are corroborated which invoke long-range Coulomb repul- 
sion as alternative explanation to understand the spectral dispersions of TTF-TCNQ quantitatively 
within an extended Hubbard model. 

PACS numbers: 71.20.Rv,61.10.Ht,71.10.Pm 



I. INTRODUCTION 

One-dimensional (ID) electron systems represent a 
valuable test case for the understanding of unusual quan- 
tum many-body states. On theoretical grounds the en- 
hanced correlations in ID lead to a breakdown of the 
paradigmatic Fermi liquid picture involving new generic 
low-energy excitations. These manifest themselves in a 
dynamic decoupling of charge and spin degrees of free- 
dom commonly referred to as spin-charge separation. 1 
While the low-energy physics is captured e. g. within the 
so-called Tomonaga-Luttinger model, which assumes lin- 
ear band dispersions and therefore does not imply any in- 
trinsic energy scale, the simplest model to explore generic 
ID physics on the energy scale of the conduction band 
width is the ID single-band Hubbard model (HM)^ The 
organic quasi-lD conductor TTF-TCNQ, where stacks 
of either the TCNQ or TTF planar molecules are effec- 
tively doped and become conducting by charge transfer 
from TTF to TCNQ, was the first material where angle- 
resolved photoemission spectroscopy (ARPES) was able 
to positively identify such ID single-band HM spectral 
features as spinon and holon branches and the shadow 
bandit However, in order to achieve a quantitative de- 
scription of the experimental data a renormalization of 
the hopping integral t at the surface by a factor of « 2 
with respect to the bulk value from density-functional 
theory or estimates from bulk-sensitive measurements 
had to be assumed. At that time it was argued that due 
to the different Madelung potential at the surface the top- 
most molecular layer may be relaxed with tilting angles of 
the molecules against the ID axis being different from the 
bulk. An alternative proposition has recently been elab- 
orated on,— namely the necessity to incorporate the next- 
nearest neighbor Coulomb repulsion V in a HM descrip- 
tion in addition to the onsite Coulomb energy U& This 



reduces the effective energy for two electrons to pass each 
other to U — V and hence enhances the effective hopping 
integral t e ff ps tU/(U — V) without the need to assume a 
larger t. Indeed, it was shown by density-functional the- 
ory that V « U/2 within a stack£ However, these results 
are based on the bulk crystal structure and an experimen- 
tal proof that no surface reconstruction is at work here 
is still lacking. There are few experimental techniques 
which in principle could give access to the molecular sur- 
face tilt such as quantitative low-energy electron diffrac- 
tion (LEED) IV measurements or scanning tunneling mi- 
croscopy (STM). Nevertheless these techniques if appli- 
cable at all provide only indirect information and need 
theoretical simulations to interpret the experimental re- 
sults. Here we employ a method widely used in surface 
physics to determine orientations of adsorbed molecules 
and apply it to the surface of TTF-TCNQ, namely angu- 
lar dependent near-edge x-ray absorption fine structure 
(NEXAFS) spectroscopy^ 



II. EXPERIMENTAL AND TECHNICAL 
DETAILS 



The NEXAFS measurements were performed at the 
PM-3 bcamlinc of BESSY (Berlin, Germany) using its 
MUSTANG endstation. The energy resolution amounted 
to 35 meV at the N K edge and was even lower according 
to the characteristics of a plane-grating monochromator 
at the C K- and S L-edges. Clean sample surfaces were 
exposed through in situ cleavage of TTF-TCNQ single 
crystals at low temperatures. To reduce radiation dam- 
age all spectra were recorded at 60 K, i. e. above the 
charge density wave (CDW) phase transition tempera- 
ture in the metallic phase. In order to achieve the re- 
quired surface sensitivity the absorption signal was mon- 
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itored in the Auger electron yield mode using a SPECS 
PHOIBOS 150 analyzer set to kinetic energies of 151 eV, 
268 eV, and 381 cV for S-, C-, and N-specific Auger lines, 
respectively. At these kinetic energies the electron mean 
free path is typically less than 10A£ which compares to 
the thickness of a single molecular layer in TTF-TCNQ 
of c/2 = 9.23 A where c is the lattice constant perpen- 
dicular to the surface (see below). Hence the absorption 
signal indeed stems essentially from the topmost molec- 
ular layer. All absorption spectra were corrected for the 
time and energy dependence of the incoming photon flux 
through division by the ring current and the drain cur- 
rent of a freshly sputtered Pt foil, respectively. The data 
were normalized in an energy region about 75 eV above 
the absorption thresholds, where the final states are ba- 
sically free-electron-likc and hence isotropic. Great care 
was exercised to avoid sample degradation upon irradia- 
tion. The irradiation exposure time of one and the same 
sample was about 40min before significant changes in 
the photoemission spectra of the S 2p core level were 
detected which turned out to monitor most sensitively 
photon-induced damages. 

In a molecular orbital picture NEXAFS involves the 
resonant x-ray excitation of a K- or L-shell electron into 
unoccupied low-lying molecular states governed by dipole 
selection rules. These project out the s-, p-, or d-like com- 
ponents of the antibonding a* or tt* final states at the 
excited atom which makes NEXAFS site- and element- 
specific. The idea to determine the molecular orientation 
of molecules on or at surfaces is based on the polariza- 
tion dependence of the involved dipole matrix elements. 
Maximum absorption is obtained when the electric field 
vector E of the linearly polarized x-rays points along a 
vector-like orbital or lies within a plane-like orbital^ The 
absorption intensity then scales as cos 2 S where 6 denotes 
the angle between E and the direction of a vector-like or- 
bital or the plane of a plane-like orbital. Detailed expres- 
sions for the angular dependence of the X-shell absorp- 
tion signal are given in Ref. considering both the polar 
and azimuthal orientation of molecular orbitals with re- 
spect to the sample surface normal and in addition the 
degree of polarization. We note that in our case the ex- 
pressions are simplified and become applicable also for 
i-shcll absorption since the azimuthal orientation of the 
TTF and TCNQ molecules can be assumed to be fixed 
and the variation of the x-ray light incidence angle is 
restricted to the plane spanned by the sample surface 
normal and the ID axis. If one moreover takes into ac- 
count that the molecules from adjacent stacks of the same 
kind are tilted in opposite directions one finally arrives 
at the following expressions for the angular dependence 
of vector and plane like orbitals, respectively: 

I v = C ■ (cos 2 9 ■ cos 2 a + sin 2 9 ■ sin 2 a) (1) 

I p = C' ■ (P ■ (1 - cos 2 9 ■ cos 2 7 - 

-sin 2 6>-sin 2 7) + (l-P)) (2) 

Here 9 is the angle between the direction of the incident 
x-rays and the sample surface, a is the angle between 



the sample surface normal and the vector pointing in 
the direction of a vector-like orbital, 7 denotes the angle 
between the sample surface normal and the normal vector 
of a plane-like orbital, P is the degree of polarization, 
and C and C are normalization constants £ Note that 
a and 7 are related to the corresponding tilt angles of 
the molecules through 90° — a and 90° — 7, respectively. 
The tilt angles can be determined by this method to an 
accuracy of about ±10°. 

In order to study the effect of surface relaxation on 
the hopping integral t band calculations were performed 
within the standard density-functional theory (DFT) us- 
ing the generalized gradient approximation (GGA)J£ We 
employed the self-consistent full-potential linearized aug- 
mented plane wave (LAPW) method as implemented in 
the WIEN2k code^i 



III. STRUCTURAL AND ELECTRONIC 
PROPERTIES OF TTF-TCNQ 

TTF-TCNQ (tetrathiafulvalene-tetracyanoquinodime- 
thane) crystallizes in a monoclinic structure (see Fig. [1]), 
space group P2\/c, with lattice parameters a = 12.298 A, 
b = 3.819 A, c = 18.468 A, and (3 = 104.46°^ It 
consists of stacks of either TTF or TCNQ molecules 
which run along the crystallographic b axis. Quasi- 
onc-dimcnsional bands are formed through overlapping 
7r type molecular orbitals which extend over the en- 
tire molecules. Maximum hopping probability along the 
stacks and hence gain in kinetic energy is accomplished 
by a tilt of the essentially planar and rigid molecules 
around the a axis by 24.5° and 34.0° for the TTF and 
TCNQ molecules, respectively. The tilt angle is oppo- 
sitely oriented between adjacent stacks. Mctallicity arises 
through charge transfer of about 0.6e per TTF to each 
TCNQ molecule. A detailed description of the strongly 
anisotropic band structure as obtained by GGA calcu- 
lations can be found elsewhere^ Here we focus on the 
TTF- and TCNQ-derived conduction bands just below 
the Fermi energy Ep along b, i. e. along TZ in the 
1st Brillouin zone (black lines in Fig. ^ . Around T 
a TCNQ-derived band doublet (split by hybridization 
of the electronic states of the two stacks per unit cell) 
can be identified while the corresponding almost degen- 
erate TTF-derived doublet is found towards the Brillouin 
zone boundary. What is important in our context is the 
fact that the experimental ARPES dispersions indicate 
a bandwidth larger by a factor of « 2 with respect to 
band theory^il This discrepancy might be reconciled by 
the assumption of a relaxed topmost molecular layer. 
At the surface the Madelung potential and the polar- 
ization screening are different from the bulk which may 
allow for different molecule tilt angles to minimize en- 
ergy. Larger tilt angles lead to a smaller intermolecu- 
lar distance d and hence an increased transfer integral 
which is counteracted by a reduced area overlap through 
a larger lateral offset I of adjacent molecules in a stack. 
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FIG. 1: (Color online) (a) Central projection view of the crys- 
tal structure of TTF-TCNQ along the b axis, (b) Parallel 
projection view along a with a hypothetical change in the 
molecule tilt angles at the surface. For details see text. 



Assuming explicitly a 1 /<i 3 -depcndenc o 13 ' 14 of the trans- 
fer integral this simple consideration leads to estimates 
for renormalized tilt angles as required to reproduce the 
experimental band widths, viz. 50° and 48.3° for TCNQ- 
and TTF-molecules, respectively. In order to explore this 
idea of a possible surface relaxation further we have per- 
formed new DFT calculations for a hypothetical bulk 
structure based on these tilt angles. While these in- 
deed almost replicate the observed larger TTF-derived 
conduction band width the occupied part of the TCNQ 
band around T comes out to have hardly any dispersion. 
Probably this is due to the complex nodal structure of 
the TCNQ-LUMO 5 perpendicular to the molecular plane 
which was neglected in the above qualitative considera- 
tion but can modulate the otherwise slowly varying func- 
tional behavior of \t\ in an oscillating manner causing one 
or several minimal 1 ^ In any case, one can induce from the 
above that (i) a molecular surface relaxation with d and 
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FIG. 2: (Color online) Bulk GGA band structure of TTF- 
TCNQ along the ID direction for the real (black solid lines) 
and a hypothetical (red solid lines) band structure based upon 
enlarged molecule tilt angles. For details see text. 



I values in a realistic parameter regime can indeed be 
responsible for the enhanced experimental band widths 
and (ii) that a sizeable increase in the tilt angles by ~ 20° 
would be necessary. We add that it was briefly noted 
in the context of recent GGA calculations that a struc- 
tural optimization for the outer molecular layer keeping 
the others fixed found changes in the molecule inclination 
angles though with minor consequences for the electronic 
dispersions along &J£ 



IV. ANGULAR-DEPENDENT NEXAFS 
SPECTRA 

The site and element sensitivity of NEXAFS enables 
the determination of the TTF and TCNQ surface tilt 
angles separately recording the S L and N K absorption 
edges, respectively. On the other hand the angular de- 
pendence of the C K edge spectra contains information 
about the orientation of both molecules provided that 
the numerous resonances can be disentangled. In order 
to accurately extract the angular dependence of a certain 
resonance feature the NEXAFS spectra have to be ana- 
lyzed using a suitable fitting scheme. We followed the ap- 
proach of Outka and Stohr— and fitted the spectra after 
a proper background subtraction with a minimal num- 
ber of Gaussian- and Lorcntzian-like peaks. Absorption 
into the free-electron continuum was accounted for by 
an error function step where appropriate, whose position 
and width reflects the binding energy and the relative 
chemical shifts due to the various inequivalent atoms in 
the structure. No "hidden" peaks were included in the 
fitting at a certain angle unless there was a definitive in- 
dication for them in the spectra at neighboring angles 
and a remnant still seen in the spectrum to be fitted. 
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FIG. 3: (Color online) Left panel: N K edge NEXAFS spectra 
of TTF-TCNQ for various light incidence angles 9. In the 
lower part the normal incidence (NI) spectrum is reproduced 
together with the corresponding fit analysis from which the 
integral peak intensities are obtained. Right panel: Integrated 
peak areas for the peaks as in the fit of the left panel as a 
function of light incidence angle (dotted curve). The solid 
lines mark fits according to the model as described in the 
text. 



Ideally, during the fitting of a complete angle series the 
parameters for all peaks (including the continuum step) 
are determined once and for all for a particular spectrum 
and then held fixed except for the ones describing the 
relative integral weight of each component. In practice, 
this cannot be sustained all over, especially for the C K 
edge where 9 crystallographically different C atoms con- 
tribute and a unique identification of each resonance for 
all angles becomes difficult. Moreover, a small relaxation 
(within a few percent) of all parameters was permitted 
in order to achieve better fits. For peak assignment we 
relied on the comprehensive NEXAFS study of Fraxedas 
et al. on highly textured TTF-TCNQ films.-^ Note that 
throughout this paper a* and 7r* resonances correspond 
to plane- and vector-like molecular orbitals within the 
classification scheme of Ref. 0. Employing Eq. [2] to fit 
the angular intensity variation of a resonance excitation 
into a plane-like orbital we assumed a polarization degree 
P fH 1 

In the left panel of Fig. [3] N K absorption spectra are 
displayed for various light incidence angles 9. On unbi- 
ased inspection four resonances can be discerned which 
are numbered and denoted according to their symmetry. 
Their intensity variations with different light incidence 
angles are clearly seen. The N K edge spectra represent 
the "clean" case, i. e. all spectra can be fitted with one 
single set of underlying components. An example of such 
a fit analysis with the total fit curve and its decomposi- 
tion into the various constituents is depicted in the lower 
part of the panel for 9 = 90°. The integral intensities of 
all four peaks are plotted in the right panel of Fig. [3] as 
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FIG. 4: (Color online) Left panel: S L edge NEXAFS spec- 
tra of TTF-TCNQ for various light incidence angles 9. Right 
panel: Integrated area for the pronounced a* resonance peak 
at 164.8 eV (dotted curve). The solid line depicts a fit accord- 
ing to the model as described in the text. 



a function of 9 together with regression curves according 
to Eqs. [T] and[2j As tilt angles for the TCNQ molecules 
one thus obtains 33.2° (34.0°, 42.8°, 28.6°) for peak 1 
(2, 3, 4). These values are in fair agreement with the 
number 34.0° from the bulk. The somewhat higher value 
derived from the intensity variation of peak 3 probably 
reflects the insufficient approximation of the underlying 
spectral weight by just one peak as can be seen from its 
relative broad Lorcntzian shape. Actually, most of the 
resonances which are denoted here comprise more than 
one molecular orbital with the same symmetries.— How- 
ever, this shortcoming has to be accepted following the 
directive to use a minimal set of peaks, for all of which 
clear indication is to be seen in at least one spectrum. 

Figure [4] shows S L edge NEXAFS spectra which are 
specific for TTF. As before distinct intensity variations 
as a function of 9 are observed for at least four resonances 
which again are marked according to their symmetry in 
the plot. A detailed fit analysis is hampered here for 
two reasons. First the signal-to-noise ratio is much worse 
in the Auger yield detection mode as compared to the 
N K edge which however cannot be easily compensated 
by longer measuring times due to the comparatively rapid 
radiation induced sample degradation. Second the over- 
all spectral shape is much broader and the individual 
resonance peaks are much less well separated than in the 
case of the N K edge spectra except for the a* resonance 
at 164.8 eV. Therefore we only fitted this single peak for 
the evaluation of the TTF tilt angle. From the fit of the 
angle-dependent integral intensity which is displayed in 
the right panel of Fig.[l]a value of 15.7° is obtained again 
in reasonable agreement with the bulk value of 24.5° and 
at variance with the expectation in case of a sizeable sur- 
face relaxation. 

The results obtained so far for the TTF and TCNQ 



5 




■ i 1 i ■ 1 i 1 ■ 1 ■ ■ 1 ■ 1 ■ ■ 1 1 ■ ■ ■ 1 ■ ■ ■ ■ 1 ■ ■ " I i I i I i I i I i I i I i I i I i L 

283 284 285 286 287 288 10 20 30 40 50 60 70 80 90 
Photon energy (eV) Angle (deg.) 



FIG. 5: (Color online) Left panel: C K edge NEXAFS spectra 
of TTF-TCNQ for various light incidence angles 9. In the 
lower part the spectrum for 9 — 20° is reproduced together 
with the corresponding fit analysis from which the integral 
peak intensities were obtained. Right panel: Integrated peak 
area for the peak marked in orange in the left panel as a 
function of light incidence angle (dotted curve). The solid 
line depicts a fit according to the model as described in the 
text. 



tilt angles might be corroborated from an analysis of the 
C K spectra if one succeeds to trace the intrinsic angular 
dependence of a resonance located either on the TTF or 
TCNQ molecules. In Fig. [5] C K spectra taken at dif- 
ferent light incidence angles are plotted. By comparison 
of C K edge spectra for TTF, TCNQ, and TTF-TCNQ 
shown in Ref. [l6| we can grossly label few structures in 
our data according to their symmetry and origin from 
either TTF or TCNQ. From the fit analysis — a rep- 
resentative regression curve is plotted in the lower part 
of the left panel in Fig. [5] — it turns out that a unique 
tracking is only possible for the most pronounced reso- 
nance structure at 285.6 eV which is identified as a it* 
resonance on TTF. In the right panel the angular vari- 
ation of its integral intensity is shown together with the 
corresponding fit curve giving a tilt angle of 27.8°. This 
again reflects the bulk value of 24.5° in agreement with 
the analysis of the S L edge spectra instead of a strongly 
increased tilt for a relaxed surface. 



V. DISCUSSION 

In the search of experimental realizations of generi- 
cally ID systems many claims have been made on the 
observation of certain isolated aspects such as spin-charge 
separation or the non-universal power-law suppression of 
spectral weight at the chemical potential. However, it is 
fair to say that a comprehensive theoretical and exper- 
imental picture where a sufficient number of pieces are 
put together thus exhibiting clear ID phenomenology has 



not yet been achieved for any of these systems. Probably 
for TTF-TCNQ such an understanding is farthest devel- 
oped. There not only the electronic ARPES dispersions 
of the TCNQ related spectral weight could be quanti- 
tatively fitted by the spinon and holon branches of the 
ID HM using the Bcthe Ansatz solution^ and the pseud- 
ofcrmion dynamical theorjsi^ but also the overall spec- 
tral weight distribution could be reproduced by density- 
matrix rcnormalization-group (DDMRG) calculations^ 
The underlying parameter set in both cases, viz. U = 
1.96 eV and t = 0.4 eV, however, assumes a hopping inte- 
gral t which is enhanced by a factor of ~ 2 with respect 
to the DFT calculations. To justify this assumption it 
was argued from an experimental point of view that the 
topmost surface layer of TTF-TCNQ might be relaxed^ 
On the other hand, further theoretical work has shown 
that the also observed transfer of spectral weight at kp 
over the entire conduction band width with increasing 
temperatures cannot be reconciled within this parameter 
seti&^Si The effective magnetic exchange J e // which de- 
fines a characteristic temperature for the onset of such 
an spectral-weight transfer would come out much too 
large within the ID HM. This questions the assump- 
tion of a surface relaxation and hence an enhanced value 
of t to explain the experimental data. Results by ex- 
act diagonalization^ for T = found that the inclusion 
of the nearest-neighbor Coulomb repulsion V in the ex- 
tended Hubbard model offers a possible loophole since 
it effectively increases the bandwidths of the spinon and 
holon excitations to match the experimental dispersions 
while at the same time allowing for a smaller t and 
hence a smaller J e // in better agreement with experi- 
ment. Indeed, Koch et alJ± recently have reported on 
density-functional results which find that the long-ranged 
Coulomb repulsion even to the third nearest neighboring 
molecule within a stack is still about 20 % of U. How- 
ever, these results are based on the bulk crystal structure 
of TTF-TCNQ. In this context the NEXAFS work pre- 
sented here provides the missing link which is needed to 
reconcile ARPES dispersions and T dependence with the- 
ory, viz. HM and DFT calculations. Our results indicate 
that no surface relaxation, i. e. no extra molecule tilt, 
occurs for either the TTF or the TCNQ stacks leaving 
the hopping integral t at the surface essentially unaltered 
with respect to the bulk. These findings validate the 
work by Koch et al. which in turn corroborates the theo- 
retical conclusion concerning the temperature-dependent 
redistribution of spectral weight above Jeff- It should 
not be concealed here that a theoretical description of 
the complete ARPES data is not yet accomplished. In 
particular, the TTF derived part of the data calls for a 
consistent treatment on the same footing. This should be 
able to explain e. g. the seeming absence of the spinon- 
holon splitting although the electronic states of the TTF 
stacks should be more strongly correlated than those of 
the TCNQ stacks. 
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VI. CONCLUSIONS 

From the analysis of the angular dependence of the 
N K, S L, and C K NEXAFS spectra we arrive at a con- 
sistent picture where under the restrictions imposed by 
experiment the topmost molecular layer of TTF-TCNQ 
single crystals appears to be structurally unaltered with 
respect to the bulk. This result lends additional signifi- 
cance to the approach to explain the deviations between 
experimental and DFT band dispersions purely electroni- 
cally. A refinement of the single-band ID Hubbard model 
by inclusion of a next-neighbor Coulomb repulsion V as 



proposed to explain the temperature dependence of the 
photoemission data thus seems plausible from an exper- 
imental point of view. 
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